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Abstract — The latest advancements in Science and 
Technology, have witnessed radical changes in the 
banking system. Today most of the banks adopt the net 
banking facility and most of the users are also addicted to 
this system. Accordingly, most of the transactions are now 
online based and much emphasis is therefore needed to 
ensure the security of authenticating a person and 
validating the transaction. Many models are therefore 
proposed in the literature. Most of these are an 
alphanumeric based password schemes or biometric 
schemes or graphical password based schemes. Each of 
these models is proposed by underlying an advantage. 
The alphanumeric passwords are proposed, with the 
assumption that generating the password is easy and the 
generated password is unique, the biometric password 
schemes are proposed with the assumption that tampering 
a biometric is next to impossible. Graphical passwords 
are proposed with an option to the user so that he can 
select an image of his choice and then select some points 
which is called a click pattern, which is unique to every 
user. 

Keywords — Graphical Password Authentication, 
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I. INTRODUCTION- PEARSONIAN TYPE III 
MIXTURE MODEL 

In the earlier research work, we developed and analyzed a 
model for segmentation based on the mixture Pearsonian 
Type I Distribution with K-means algorithm. These 
models are useful when the pixel intensities of the feature 
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vector in the image regions are left skewed. But in some 
image the pixel intensities of the image regions may not 
be left skewed. They may have long upper tail with a 
right skewed nature. To segment these types of images, it 
is needed to consider that the pixel intensities of the 
image regions follow a right skewed distribution. Hence, 
in this research article we develop and analyze image 
segmentation method based on mixture of Pearsonian 
Type III Distribution. The Pearsonian Type III 
Distribution is capable of portraying right skewed and 
long upper tail distributions. Here, it is assumed that the 
pixel intensities of each image region follow Type III 
Pearsonian Distribution and the pixel intensities in the 
whole image is characterized by a finate mixture of 
Pearsonian Type III Model. 

II. PEARSONIAN TYPE III MIXTURE 
DISTRIBUTION 

In any image analysis, the entire image is considered as a 
union of several image regions. In each image region the 
image data is quantified by pixel intensities. The pixel 
intensity z = f( x , y ) for a given point (pixel) (x, y) is a 
random variable, because of the fact that the brightness 
measured at a point in the image is influenced by various 
random factors like vision, lighting, moisture, 
environmental conditions etc. To model the pixel 
intensities of the animal and human image regions, it is 
assumed that the pixel intensities of the region follow a 
Pearsonian Type III Distribution (PTIIID). The 
probability density function of the pixel intensity is 


< 00,-00 < q y < 00 (Equation-1) 


Where □ is a Gamma Function 
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For different values of the parameters the various shapes 
of probability curves associated with Pearsonian Type III 
Distribution are shown in Figure-1. 



Fig.l: Frequency Curves of Pearsonian Type III 
Distribution 

The entire animal and human image is a collection of 
regions which are characterized by Pearsonian Type III 
Distribution. Here, it is assumed that the pixel intensities 
of the whole image follow a K-component mixture of 
Pearsonian Type III Distribution and its probability 
density function is of the form: 

K 

p{z)=Y J a if i {^ a i^q l ) 

i=l 

Where, K is number of regions, 0 < « ; < 1 are weights 

such that a i = 1 and is ft (-'■ I a i > *7/ ) as given in 

equation (1) a i is the weight associated with i th region 
in the whole image. 

In general the pixel intensities in the image regions are 
statistically correlated and these correlations can be 
reduced by spatial sampling (Lei T. and Sewehand W. 
(1992)) or spatial averaging (Kelly P.A.et al (1998)). 
After reduction of correlation, the pixels are considered to 
be uncorrelated and independent. The mean pixel 

K 

intensity of the whole image is E(Z)=Z««« 

i=l 

m. INDENTATIONS AND EQUATIONS 

In this section, it is derived estimates of the model 
parameters using Expectation Maximization (EM) 
algorithm. The likelihood function of the observations 
given by 


Zi, Z 2 ,....,Z N drawn from an image with probability 


density function 
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(Equation-2) 

The beginning step of the Expectation Maximization 
algorithm requires the initialization of two parameters 

( a i e Ji ; i=l,2,. . .k) and component weights 

(cr, , i — 1,2,... ,K) f rom t| le observed values. The idea 
of EM algorithm is to iteratively calculate Maximum 
likelihood estimates of unknown parameters 
6 = (a i q i a i \i = \,2,...,K) 


E-STEP : 

In the expectation (E) step, the expectation value 
of log with respect to the initial parameter vector 

e m 

Q(6>;0 m )=E 0W [logL(6>)/z] 

Given the initial parameters 0 " r) , one can compute the 
density of pixel intensity z s as 

i=l s=l 

K K 

This implies logL(#) = £log £a i /,(z s ,0 ,1, 

i=l L i=l 

(Equation-3) 

The conditional probability of any observation z s , 
belongs to any region K is 


t k (z s ,0 a A = 


\ cx k m f k tz s ,0 m ~\ 


oc k m f k tz„0 m 

Ptz s ,0 m ) 


^rz k w f k tz s ,0"A 
_ 1=1 


The conditional of the log likelihood function of the 
sample is 

Q(d-d m )=E 0<a \ogLW)ll] 
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Following the heuristic arguments of Jeff A. Bilmes 
(1997) we have 

Q(0-0 m )= ||((,( z „9 |1| )Xk)g fjz s ,0 m )+ log «/'>)) 

(Equation-4) 

But we have 


\ (qfii+l) 
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and 


For obtaining the estimation of the model parameters one 

has to maximize q{o-, o m ) such that 2: = 1 . This 

can be solved by applying the standard solution method 
for constrained maximum by constructing the first order 
Lagrange type function, 

,* = r^Clog » - 1 - J^l — '5 T.cc, w J 

(Equation-5) 

Where, is Lagrange multipliers combining the 
constraint with the log likelihood function to be 


q(&; 6» ll) ) = jr (t ‘ (A, , )((log f(z s , o'" ) + log a,- (n )) maximized 
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Summing both sides over all observations, we get 
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For updating the parameter ; i=l,2,. . .k we consider the derivative of 
with respect to a i and equate it to zero. 

Q(0-0 m )=E[(logL(0 m )\ 


Q(010 W ) 


We have 


Therefore 


implies is[(logL(# (,) )]= 0 


6a, 


q(6-0^)= 0 
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For updating the parameter q,, i=l,2,....,K we consider the derivative of Q 
With respect to <7, and equate it to zero 
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IV. K-MEANS CLUSTERTING ALGORITHM 

The K-means algorithm is one of the simplest clustering 
technique for which the objective is to find the partition 
of the data which minimizes the squared error or the sum 
of squared distances between all points and their 
respective cluster centers (Rose H. Turi, (2001)). K- 
means algorithm uses an iterative procedure that 
minimizes the sum of distances from each object to its 
cluster centroid, over all clusters. This evaluational 
procedure consists of the following steps. 

1. Randomly choose K data points from the whole 


data set as initial clusters. These data points 
represent initial cluster centroids. 

2. Calculate Euclidean distance of each data point 
from each cluster center and assign the data points 
to its nearest cluster center. 

3. Calculate new cluster center so that squared error 
distance of each cluster should be minimum. 

4. Repeat step II and III until clustering centers do 
not change. 

5. Stop the process. 

In the above algorithm, the cluster centers are only 
updated once all points have been allocated to their 
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closed cluster center. The advantage of K-Means 
algorithm is that it is a very simple method, and it is 
based on intuition about the nature of a cluster, which is 
that the within cluster error should be as small as 
possible. The disadvantage of this method is that the 
number of clusters must be supplied as a parameter, 
leading to the user having to decide what the best number 
of clusters for the image is (Rose H. Turi, (2001)). 
Success of K-means algorithm depends on the parameter 
K, number of clusters in image. 

After determining the final value of K (number of 
regions), we obtain the initial estimates au q h and «, of 
for the ith region using the segmented region pixel 
intensities with Pearsonian Type III Distribution. The 
initial estimate a,, a,=l/K, where is taken as i=l,2,...,K. 
The parameters a, and q, are estimated by the method of 
moment’s pi as first moment pi and its three central 
moment’s P 2 P 3 and p 4 . 

V. CONCLUSION 

The number of image regions is estimated by utilizing the 
image histogram. The model parameters are estimated by 
deriving the updated equations of EM algorithm. The 
initialization of the model parameters is done through K- 
means algorithm and moment method of estimation. The 
segmentation algorithm is developed through maximizing 
the component likelihood under Bayes framework applied 
in MIR Flickr . 
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